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1. Introduction 

Low-energy QCD can be characterized by four nonperturbative quantities, of which the gluon 
condensate ({gF^ v ) 2 } and the vacuum correlation length X (that is, the distance at which the two- 
point, gauge-invariant, correlation function of gluonic field strengths exponentially falls off) are 
related to confinement (so that the string tension corresponding to the two static sources in the fun- 
damental representation is o °< X 2 ((gF£ v ) 2 ) [|l|]), while the quark condensate (y/'y), together with 
the constituent quark mass m, characterizes spontaneous breaking of chiral symmetry. The quark 
condensate is expressible in terms of gluonic degrees of freedom, (yy) = —d(r[A a ^,m])/dm, 
where (r[A^,ra]), the one-loop effective action, can be represented as a world-line integral over the 
closed quark's trajectories {z^(t)} and their anticommuting counterparts {y^i(t)} [0], 



(r[A a u , m])=-2Nf n- e ->" 2 > [ @ z f ®y ■ e SZd<z*l+in%) 
F Jo s Jp Ja 



X 



tr £P exp 



ig J' W a (A£z M - Wv F£ v ) 



N c }. (1.1) 



In ( |l.l| ), Nf is the number of light-quark flavors, s is the Schwinger proper time "needed" for the 
quark to orbit its Euclidean trajectory, F^ v = e^A" — d v A^ + gf ahc A^A^ is the Yang-Mills field- 
strength tensor, and T a 's are the generators of the group SU(A^ C ) in the fundamental representation, 
obeying the commutation relation [T a ,T b ] = if" bc T c . Notice that, since the quark condensation is 
generally argued to occur due to the gauge fields, the free part of the effective action in Eq. ( |iH| ) 
has been subtracted, so that (r[0,m]) = 0. Furthermore, P and A stand there for the periodic (J p = 
J Zfl ( s )= ZM (o)) an d the antiperiodic (J A = fy (j) = _y (o)) boundary conditions, which are imposed, 
respectively, on the trajectories z^(t) and their Grassmannian counterparts Yfi( T ) describing y- 
matrices ordered along the trajectory. The trajectories obey the condition Jq dzz^(T) = 0, meaning 
that the center of each trajectory is the origin. That is, the factor of volume associated with the 
translation of a trajectory as a whole is divided out, and the vector-function z^t) describes only 
the shape of a closed trajectory, not its position in space. Finally, throughout this talk, we mean 



by the quark mass the minimal value of the mass parameter m, entering Eq. (|1.1|), which renders a 
finite (\jf\j/) — see Ref. [||]. 

From the mathematical viewpoint, the effective action ( |1.1[ ) represents an integral over closed 
quark trajectories, with the minimal surfaces bounded by those trajectories appearing as arguments 
of the Wilson loops. That is, 

(r[A£,m]> = -2N f £ '^-e-™ 2 * J @ v ^- fo d <Ul+^w) x 

x {exp [-2[d^ WvJ -^ (W[ Z ,])-N C ] , (1.2) 
with the Wilson loop given by 

(W%]> = ^tr^ exp (ig J* d-uT^z,^ \ (1.3) 

and 8 / 8s ^ v being the area derivative operator, which allows us to recover the spin term ~ y v F^ v 
in Eq. (|L~I|). Equation (1.2) allows us to reduce the gauge-field dependence of Eq. ( |1.1[ ) to that of a 
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Wilson loop. The Wilson loop is unambiguously denned by the minimal-area surface bounded by 
the contour 

World-line integrals of this type were first calculated by imposing for the minimal surface a 
specific parametrization, which in 3D corresponds to a rotating rod of a variable length [Q]. We 
notice that similar parametrizations of minimal surfaces spanned by straight-line strings (which 
interconnect quarks or gluons in the one-loop diagrams) are widely used in the existing litera- 
ture [|J. With the 4D parametrization of this kind, we were able to obtain, within the effective- 
action formalism, a realistic lower bound for the constituent quark mass: m = 460MeV for (Yy) — 
— (250MeV) 3 (cf. Ref. [^]). For bookkeeping purposes, an explicit derivation of this parametriza- 
tion from the Nambu-Goto string action will be provided in Section II of this talk. 

It happens that for sufficiently large m's, the mean size of the trajectory is smaller than the 
vacuum correlation length A, so that the nonperturbative Yang-Mills fields inside the trajectory 
can be treated as constant, leading to the area-squared law for the Wilson loop [§]. By using 
the world-line representation (|lj]), the area-squared law can be shown (see Ref. ^ for details) 
to yield the known heavy-quark condensate of QCD sum rules, (YY)heavy x — {(s^uv) 2 ) / m > i n 
which case m becomes the constituent mass of a heavy quark. For lighter quarks, the mean size of 
their Euclidean trajectories exceeds A, so that one can expect the area-squared law to be morphed 
into an area law. However, unlike static color sources transforming according to the fundamental 
representation of SU(iV c ), which yield an area law with a constant string tension a, light quarks are 
subject to a zigzag-type motion. This type of motion can only be reconciled with the area law for 
some effective scale-dependent string tension cr(j) oc l / s , so as to obtain a nonvanishing (YY) t@]- 
We notice that the fractalization of quark trajectories, which takes place upon the deviation from the 
heavy-quark limit, has been studied in Ref. [Q] in terms of the velocity-dependent quark-antiquark 
potentials. 

It is then clear from the above discussion that the minimal value for the mass parameter m 
and the mean size of the quark trajectory are not independent from each other so that, in general, 
(r[A" m]) is a functional of the minimal area. To this we have to add that the quark condensate 
characterizes the QCD vacuum, and, therefore, should remain invariant under the variations of 
the minimal area. The final ingredient should come from energy conservation: the total energy 
of a given excited quark-antiquark system should be equal to the sum of the quark masses and 
the energy stored in the quark-antiquark string. Then, the purpose of our analysis is to study the 
effective action as a functional of the minimal area for quark-antiquark excited systems, that is, 
to evaluate the corresponding variation 8m such that (YY) remains constant for a given series of 
quark-antiquark radial excitations. At this point, we need a physical input on how to quantify these 
excitations. A natural choice will be to link these excitations to the daughter Regge trajectories [gj. 
In this work, we adopt this point of view. 

The talk is organized as follows. In the next Section, using as an input the excitation energy 
corresponding to the n-th daughter Regge trajectory, we introduce an Ansatz for the scaling factor 
Z„, which describes an increase in the area of the excited-string world sheet. In this way, we 
consider only the radial excitation modes of the quark-antiquark pair, which correspond to the 
"breathing modes" of the string world sheet. Consideration of angular excitations, which would 
correspond to the disclinations of the world sheet, lies outside the scope of our present analysis. 
Then, we proceed to a selfconsistent determination of the "critical index" y, which defines that 
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Ansatz. With this knowledge at hand, we calculate the constituent quark-mass correction 8m n as 
a function of the radial-excitation quantum number n. In Section III, using the large-71 asymptotes 
of the formulae obtained for Z n and 8m n , we show that the primary contribution to the excitation 
energy of the quark-antiquark pair stems from the constituent quark mass m n , and not from the area 
increase. We also obtain the lowest (that is, n = 1) correction to the constituent quark mass, 8m\, 
which turns out to be about 26 MeV. Also in Section III, we present some concluding remarks and 
an outlook. 



2. A correction to the constituent quark mass 

In this Section, we derive a general expression for the correction to the constituent quark mass, 
coming from radial excitations of the quark-antiquark string sweeping the surface of the Wilson 
loop. To this end, we use for the eigenenergies of radial excitations of the quark-antiquark pair the 
Regge formula [||] : 

E„ = y/no(4n + 3). (2.1) 

Here, n is the quantum number of a radial excitation, and G ~ (440MeV) 2 stands for the string 
tension in the fundamental representation of the group SU(3). 

The energy gap E„ — Eq can be filled in by both deformations of the quark-antiquark string 
and/or by increasing in the quark constituent masses. That is, 

E n -E = o ■ (L n -L ) + 2(m n -m ), (2.2) 

where L„'s are the eigenvalues of the length of the string. Notice that m n turns out to be an implicit 
function of L n . We furthermore denote by 2R n the diameter of the semiclassical Euclidean trajec- 
tory performed by the quark in the n-th excited state of the system. The value of 2R n cannot exceed 
some 2/?„ max , at which string-breaking occurs. The string elongation is given by the ratio L u /Lq, 
where in the (n = 0) state we must have 2Rq = Lq. 

Since we will be calculating the constituent quark masses m„ in the units of ^/o, it is conve- 
nient to introduce a dimensionless function 

f n = m n /y/%a. (2.3) 

The sought correction to the constituent quark mass can be obtained from the effective action ( p~T| ) 
if one uses there for the surface, entering the Wilson loop, the world sheet of the excited quark- 
antiquark string. To this end, let us start by defining the scaling factor Z„ which describes an 
enlargement of the world-sheet area over its value in the (n = 0) state. We have 



— = 1 +S n , where S„ = — 



v^+3-\/3-2-5/ n 



(2.4) 



and 8f„ = f n — fo. Accordingly, the scaling factor, which describes an increase of the area of the 
string world sheet in the n-th excited state, reads 

Zn=(^) =(1+S„) 2 . (2-5) 



L 
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We further notice that the role of a global characteristic of the quark trajectory can be played 
either by its semiclassical radius R n , or by the proper time s during which the trajectory is orbited 
by the quark. Therefore, these two quantities are related to each other through a scaling relation of 
the form 

( s V 

Rn —Rn,max ' I I • (2.6) 

V^max / 

Here s max is the proper time necessary for the quark to orbit a trajectory of the maximum diameter 
2/?« max equal to the string-breaking distance. The actual value of the "critical index" y will be 



selfconsistently determined below. Thus, we get for the correcting term S n in Eq. ( 2.5 ) the following 
expression: 



5„ = 5 n (,) = % where §, ee s ™ [v4^+3 - V3 - 2 • 8f n 



(2.7) 



We proceed now to the discussion of parametrizations for the minimal area of the string world 
sheet and for the Wilson loop. For the minimal area of the unexcited-string world sheet we use the 
following parametrization: 

1 f s . 

54D = ^7f7 dT \ e nvX P ZxZp\- (2-8) 

It represents a four-dimensional generalization of Ssu = ^ J s dl\z x ±\ (cf. Ref. [|j]), which is the 
area-functional of a surface swept out by a rotating rod of a variable length. We notice that S4U 
stems directly from the usual formula for the area (corresponding to the Nambu-Goto string action) 
upon the parametrization of the surface by the vector-function Wm(^i,Q = £2 ■z^{^\/o), where 



£1 = ctt and £2 6 [0, 1]. Indeed, the usual formula for the area reads = J s dC,\ J dtjz \/detg, 
where gab = d a w^ ■ dyw^ is the induced-metric tensor, and each of the indices a and b takes the 
values 1 and 2. Using the above parametrization for w^Cij £z)> one can then readily prove the 
following equality: 

S 4D = d = 2 / d *\jz 2 iii 2 v-{znZii) 2 - 

Next, we follow Ref. [fl for what concerns the parametrization of the Wilson loop ( |1.3| ), which can 
be written in the form 

{W[Z ^ ]) = 2»~T(a) ' (*M a - K «(*M- (2-9) 

Here, L^ v = £^ v ip Jo dlzxip is the integrated surface element, whose absolute value is implied in 

/ \ 1/2 

the sense that |£^ v | = (2 £ I?, v . Furthermore, T(x) and K a (x) in Eq. ( |2.9[ ) stand, respectively, 

for Gamma- and MacDonald functions, a > 1 is some parameter, and a = d(s) stands for the 
effective string tension of dynamical quarks. As it was shown in Ref. [^J], Eq. provides an 
interpolation between the area law for large loops and the area-squared law [|]] for small loops. 
This statement is illustrated by Figs. 1 and 2, which compare the combined parametrization (^9|), 
for a = 1.9, with the area and the area-squared laws. For this illustration, we choose a circular 
contour with the minimal area S, set N c = 3, and use the relation [^] ((gF^ v ) 2 ) = 1 -^ jj, where 
X = 1.72GCV- 1 |g] for the case of QCD with dynamical quarks considered here. The chosen 
value of a = 1.9 has been shown in Ref. [0] to provide the best analytic approximation of the 
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0.5 1 1.5 2 2.5 3 3.5 4 



Figure 1: 2 a-i' r ( a ) ■ (eS) a ■ K a (aS) at a = 1 .9, e~ aS , e 5 ~ . The value aS = 4 coiTesponds to the 

radius of the contour equal to 0.51 fm. 



'Combined' parametrization ^™ 

Area law 

'Area-squared' law 




Figure 2: In 



(aS) a -K a (aS) 



at a = 1 .9, aS, <( ^ )2) S 2 . The value aS = 5.5 corresponds 



_2 a -T(a) 

to the radius of the contour equal to 0.6 fm. The gap between — ln[- • •] and aS, being ~ ln(<r5') at large 
distances, is irrelevant for the static potential. 



area-squared law by the combined parametrization (L9) at small distances. Figure 1 illustrates the 
efficiency of this approximation. 

We are now in a position to calculate (T[A" , m n ]), which is defined by Eq. with m replaced 
by m n . For this purpose, we multiply the infinitesimal surface element dxz\Z p in Z^ v above by the 
scaling factor Z„. Following Ref. [||], we arrive at the expression 
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where 



7(Z„,^ V )= [ I 

.IP JA 



(2.11) 



and & 



The world-line integral (2.11) describes an infinite sum of one-loop quark diagrams, each 
having its own number of external lines of the auxiliary gauge field B^ v . Similarly to what was 
done in Ref. [Qj, we retain in this diagrammatic expansion only the two leading terms - the one 
corresponding to a free quark, which cancels out in Eq. (LI), and the one corresponding to a 
diagram with two external lines of the gauge field. The latter yields the following expression: 



where B 



1/2 



1 



(4ns) 7 



' : (sBZ n ) 2 + 0((sBZ n ) 4 ) 



(2.12) 



We start our analysis with the limit «> 1, where one can approximate the factor Z„ in 
Eq. ( 2.12 ) by S„. In the same large-rc limit, it is legitimate to disregard the quartic and the higher 
terms in the expansion ( 2.12| ), provided the amplitude B is bounded from above as 

1 s 2 r- 1 



B< 



(2.13) 



In the last equality, we have used the explicit parametrization ( J2.7D of S„ in terms of s. Now, as 
it was already mentioned in Introduction, the quark condensate, {\jf\jf) = — gff (r[A^,mJ), being 
the quantity which characterizes the QCD vacuum, should remain ^-independent. For this to be 
possible, the energy of string excitations should be largely absorbed by the constituent quark masses 
m n , and it will be demonstrated below that this is what indeed happens. 

An expression for the quark condensate following from Eqs. ( 2.1Q )-( 2.13 ) reads 



(W) 



a(a + l)(a + 2)Nf 



m„ 



ds- 



dB- 



B< 



(>+&) 



a+3 ' 



(2.14) 



where from now on we set A^ c = 3. The B-integration in this formula can be performed analytically, 
to yield 

f[A n (s),a] 



(W) = -J=2- m n 



3N f 
An 2 



where 



dse 



,27-1 



2s 2 A n (s) 



(2.15) 



(2.16) 



while a somewhat complicated function f[A, a] was introduced in Ref. [0]. An important property 
of this function is that, for a > 1 of interest, the ratio f[A, a] /A, as a function of A, has a maximum 
at A <C 1, which sharpens and reaches a finite value (~ 1.18), with the further increase of a. 

In order for the quark condensate to stay finite in the small-mass limit, one should be able to 
represent 



/[A " W '°' 1 in the form (cf. Refs. [jfij g|) 



f[A n (s),a] a, 



3/2 



2^ 2 A„( 



(2.17) 
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where Oq is some parameter of dimensionality (mass) 2 unambiguously related to the phenomeno- 
logical value of (YW)- Moreover, this representation should remain valid up to the values of the 
proper time such that 



m n s max > 1. 



Equation (2.17) can equivalently be written as 

/[•An, a 



x, where x = 2(oqs) 



3/2 



(2.18) 



(2.19) 



Owing to the above-mentioned form of the function f[A,oc]/A, a solution to Eq. ( |2.19| ), which 



provides a physical decrease of a with s, reads A n ~ x e , where e — > for a > 1 of interest, and 



x < 1 (cf. Ref. [|]). Therefore, to a very good approximation, one can set in Eq. ( |2.16D A n ~ 1 . 



This allows us to obtain the actual value of the power y in the initial Ansatz (2.6). To this end, 
we notice that, due to the energy conservation in the quark-antiquark system, a variation of the 
semiclassical radius R n of the trajectory leads to the variation 8m n = d8R„ of the constituent 
quark mass. Therefore, a difference between the values of the radius R n in the n-th and the 0-th 
states, 8R n =R n —Ro, reads 



8R„ 



(2.20) 



and ( 2.2C ) to have the 



a - s 2 ?- 1 ' 
where Eq. ( 2.16| ) has been used at the final step. Now, in order for Eqs. 
same ^-dependence, y= 1/3. 

With this value of y at hand, we can now calculate the correction 8m n to the constituent quark 
mass, which is produced by the radial excitations of the quark-antiquark string. To this end, we 
first insert the value of y= 1/3 into Eq. ( |2.20| ), that leads to the following relation: (8R n ) max = 
\fl8m n • ^„ 2 max Smax- Furthermore, for ^„, max in this formula we use its expression provided by 
Eq. (O) with 7=1/3. That yields 



(8R n )r 



nV2- 



8tn n s max 



n + 



(8fn)r 



(2.21) 



Next, we use the approximation 8m n s max ~ l/(5/n n ) mm , which reflects the fact that the trajectory 
of a maximum size (and therefore requiring the maximum proper time to be orbited) is reached 
when the value of the constituent quark mass is minimal. [Notice that this approximation paral- 



lels condition fl2.18| ).] Substituting this approximation into Eq. ( [2.2 H ), we arrive at the following 
equation: 

7T\/2 



(8m n ) 



^lmJSR n ) 



n + - 



(Sfn) 



(2.22) 



In order to solve this equation, we represent it entirely in terms of (5/„) m i n . That can be done by 
virtue of the relation (5m„) mm = *f%G (8f n ) m i n , which stems from Eq. ( |2.3[ ). As a result, we obtain 
the following quadratic equation: 



(8f n )mm + b- (Sf„) 



1/2 
min 



0, where b ; 



o 3 ^Ro, ma ,(SRn 

(2^)V4 



,1/2 
i max 



(2.23) 
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A solution to this equation yields the sought correction to the constituent quark mass: 

r I i 2 

(8m n ) 



f no-{8f n ) 




The limits of this formula at large and small n's will be analyzed in the next Section. 



(2.24) 



3. The limiting cases of large and small excitations. Concluding remarks 



Ignoring for a moment the effect of string-breaking, we see that the obtained (5m„) mm , Eq. 
(2.24), vanishes in the limit of b— which corresponds to the quark-antiquark string of an infinite 



length [cf. the definition of b in Eq. ( 2.230 ] . In reality, however, the string-breaking phenomenon 



imposes an upper limit on the possible values of b. Indeed, the upper limit for both /?o,max and 
(8R n ) max is given by d & ,\,,/2, where d s ,b, is the string-breaking distance. Lattice simulations and 
analytic studies [11] suggest for this distance the value of d s .\,. ~ 1.5fm. Using also the phe- 
nomenological value of a ~ (440 MeV) 2 , we get b < 1.37. Therefore, we find from Eq. ( 2.24 ) the 
realistic asymptotic behavior of the constituent quark mass to be 



(8m n ) r 



\Jn<3n for n^> I. 



(3.1) 



Comparison of this result with the initial Eqs. (|2JJ) and ( |2.2| ) shows that, in the large-w limit, 
the leading contribution to the excitation energy E„ of the quark-antiquark pair stems from the 
constituent quark masses. Indeed, one can perform the large-w expansion of (5/„) mm given by 
Eq. ( 2.24 ), which yields (5/„) mm = y/n [l — bn~ l l A + ^(n -1 / 2 )] . Then, inserting this expansion 
into the formula for S n , Eq. (|2.7|), we obtain the leading large-w behavior 



^0, max 



■bn Y '\ 



(3.2) 



which is subdominant compared to Eq. (3.1). Recalling Eq. (2.4), we conclude that 



m 



a. mm ~ L n . 



Thus, the constituent quark mass appears as a primary ingredient of the excitation energy of the 
quark-antiquark pair in the large-rc limit, whereas the elongation of the string plays only a secondary 
role. Still, we observe an increase of S n with n, which, for sufficiently large n's, validates the 
approximation Z„ ~ S 2 used after Eq. ( 2.12 ). 



Let us now evaluate a correction to the constituent quark mass, which is associated with the 
(n = 1) excitation. We note that this excitation is developed just on top of the maximally-stretched 
unexcited-string configuration. For this reason, one can use for such an evaluation the above- 
adopted maximum values of b = 1.37 and /?o,max = 0.75 fm, and also set s = s max . Extrapolating 

1.45. The fact that this extrapolation to n = 1 
1, signals the need to introduce 
Zi, where 

' 2 This equation yields the value of k ■■ 



then Eq. (3.2) down to n = 1, we get S\ 



SiL, 

1 J — " max 

of the initial parametrically large result of Eq. (3.2) leads to S\ 
some correcting numerical factor k. It can be defined through the relation {kS\) 2 
again Z\ = (1 +S\) 2 . This equation yields the value of k = 1.69. Next, according to Eq. (3.2), the 
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multiplication of S\ by a factor of k is equivalent to the multiplication of b by such a factor. This 
observation yields the corrected value of b = 2.32. Inserting it into Eq. ( |2.24 ), we get the sought 
estimate for the [n = 1) correction to the constituent quark mass: 

(<5mi) min = 26.0MeV. 

This value looks like a reasonable additive correction to the leading result, m = 460 MeV, quoted 
in the Introduction. 

In conclusion, we notice that excitations of the quark-antiquark pairs can in general lead to 
an increase of the constituent quark mass and to an elongation of the quark-antiquark string. In 
this talk, we have shown that, for large radial excitations n, the constituent quark mass grows 
as &(n l l 2 ), while the length of the string grows only as 0(n l l A ). This result clearly means that 
the excitation energy of a quark-antiquark pair stems mostly from the increase of the constituent 
quark masses and not from string elongation. Thus, at least within the effective-action formalism, 
excited quark-antiquark bound states tend to have essentially the same size, irrespective of their 
radial excitations. This is true even if we had energy dependence for bound states different from 
that of Eq. (|Q|), because that will be just another prescription for the area scaling and hence, the 
final qualitative result would not depend on actual details on how this scaling is obtained, but just 
from the fact that for a given scaling up of the area, the mass would go like the square of the 
string elongation, whereas the string energy would dimensionally go with oL n . Finally, such size 
stiffness will preclude decay channels to become large with string elongations, because they chiefly 
measure the size of the parent hadron and that does not change appreciably. We also emphasize 
that the adopted calculational method does not rely on any specific class of the string deformations 
(such as, e.g., the normal modes). Finally, it looks natural to apply the present approach to the 



description of an interesting lattice result [12] that, for quarks in the fundamental representation, 
the deconfinement and the chiral-symmetry-restoration temperatures are nearly the same, whereas 
for quarks in the adjoint representation, the chiral-symmetry-restoration temperature exceeds the 
deconfinement one by a factor of 8. Work in this direction is currently in progress. 
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